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Abstract 

This paper is devoted to investigate the exact solutions of Bianchi 
types I and V spacetimes in the context of f(R,T) gravity pQ. For 
this purpose, we find two exact solutions in each case by using as- 
sumption of constant deceleration parameter and the variation law of 
Hubble parameter. The obtained solutions correspond to two different 
models of this universe. The physical behavior of these models is also 
discussed. 

Keywords: f(R,T) gravity, Bianchi types I and V, deceleration param- 
eter. 

PACS: 04.50.Kd, 98.80.-k, 98.80.Es. 



t 



1 Introduction 

The most popular issue in the modern day cosmology is the current expansion 
of universe. It is now evident from observational and theoretical facts that our 
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universe is in the phase of accelerated expansion [2] -[ID]- The phenomenon 
of dark energy and dark matter is another topic of discussion [11J-[18J. It 
was Einstein who first gave the concept of dark energy and introduced the 
small positive cosmological constant. But after that, he remarked it as the 
biggest mistake in his life. However, it is now thought that the cosmological 
constant may be a suitable candidate for dark energy. Another proposal to 
justify the current expansion of universe comes from modified or alternative 
theories of gravity. f(T) theory of gravity is one such example which has been 
recently developed. This theory is a generalized version of teleparallel gravity 
in which Weitzenbock connection is used instead of Levi-Civita connection. 
The interesting feature of this theory is that it may explain the current 
acceleration without involving dark energy. A considerable amount of work 
has been done in this theory so far [19] . Another interesting modified theory 
is f(R) theory of gravity. In this theory, a general function of Ricci scalar is 
used in standard Einstein-Hilbert lagrangian. Some review articles [20J can 
be helpful to understand the theory. 

Many authors have investigated f(R) gravity in different contexts. Spher- 
ically symmetric symmetric solutions are most commonly studied solutions 
due to their closeness to the nature. Multamaki and Vilja [21] explored vac- 
uum and perfect fluid solutions of spherically symmetric spacetime in metric 
version of this theory. They used the assumption of constant scalar curvature 
and found that the solutions corresponded to the already existing solutions in 
general relativity (GR). Noether symmetries have been used by Capozziello 
et al. [22] to study spherically symmetric solutions in f(R) gravity. Sim- 
ilarly many interesting results have been found using spherical symmetry 
in f(R) gravity [23]. Cylindrically symmetric vacuum and non- vacuum so- 
lutions has also been explored in this theory [21]. Sharif and Shamir [25J. 
found plane symmetric solutions. The same authors [2S] discussed the so- 
lutions of Bianchi types I and V cosmologies for vacuum and non-vacuum 
cases. We [27] calculated conserved quantities in f(R) gravity via Noether 
symmetry approach. 

In a recent paper [lj, Harko et al. proposed a new generalized theory 
known as f(R, T) gravity. In this theory, gravitational Lagrangian involves 
an arbitrary function of the scalar curvature R and the trace of the energy- 
momentum tensor T. Myrzakulov [28J discussed f(R,T) gravity in which he 
explicitly presented point like Lagrangians. The exact solutions of f(R,T) 
field equations for locally rotationally symmetric Bianchi type I spacetime 
has been reported by Adhav [2H] • Sharif and Zubair [5U] discussed the laws of 
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thermodynamics in this theory. Houndjo [31 J reconstructed f(R, T) gravity 
by taking f(R,T) = fi(R) + fz(T) and it was proved that f(R,T) gravity 
allowed transition of matter from dominated phase to an acceleration phase. 
Thus it is hoped that f(R, T) gravity may explain the resent phase of cosmic 
acceleration of our universe. This theory can be used to explore many issues 
and may provide some satisfactory results. 

In this paper, we are focussed to investigate the exact solutions of Bianchi 
types I and V spacetimes in the framework of f(R, T) gravity. The plan 
of paper is as follows: In section 2, we give some basics of f(R, T) gravity. 
Section 3 and 4 give the exact solutions for Bianchi types I and V spacetimes. 
Concluding remarks are given in the last section. 



2 Some Basics of f(R, T) Gravity 

The f(R, T) theory of gravity is the generalization or modification of GR. 
The action for this theory is given by p] 

s = J V=g (j^ffa t ) + L -) d * x > C 1 ) 

where f(R,T) is an arbitrary function of the Ricci scalar R and the trace T 
of energy momentum tensor T^ u while L m is the usual matter Lagrangian. It 
is worth mentioning that if we replace f(R, T) with f(R), we get the action 
for f(R) gravity and replacement of f(R, T) with R leads to the action of 
GR. The energy momentum tensor T^ u is defined as 



rp 2 5(^gL m ) 

v ~ sir ■ {2} 

Here we assume that the dependance of matter Lagrangian is merely on the 
metric tensor g^ u rather than its derivatives. In this case, we obtain 

T^ v = L m g^ - 2 — ^. (3) 

The f(R,T) gravity field equations are obtained by varying the action S 
Eq.flT]) with respect to the metric tensor 



f R {R,T)R^~f{R,T)g^-{V^ u -g^U)f R {R,T) = kT^-MR, T)(T^+0^), 

(4) 
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where denotes the covariant derivative and 

n _™^7 { i p rp\ df R (R,T) df R (R,T) a ^ T ^ 

Contraction of (jlj) yields 

f R (R, T)R + 3D/ fl ( J R, T) - 2f{R, T) = kT - f T (R, T)(T + 0), (5) 

where G = M M . This is an important equation because it provides a rela- 
tionship between Ricci scalar R and the trace T of energy momentum tensor. 
Using matter Lagrangian L m , the standard matter energy-momentum tensor 
is derived as 

T„ v = (p + p)u,,u v - pg^ u , (6) 

where = ^^00(1,0,0,0) is the four-velocity in co-moving coordinates and 
p and p denotes energy density and pressure of the fluid respectively. Perfect 
fluids problems involving energy density and pressure are not any easy task 
to deal with. Moreover, there does not exist any unique definition for matter 
Lagrangian. Thus we can assume the matter Lagrangian as L m = —p which 
gives 

= -V9ia> ~ 2T ilv , (7) 
and consequently the field equations (jlj) take the form 

f R (R, T)Rp,~f(R, T)g^-(V^„-g^n)f R (R, T) = KT^+f T (R, T)(T^+pg^), 

(8) 

It is mentioned here that these field equations depend on the physical nature 
of matter field. Many theoretical models corresponding to different matter 
contributions for f(R,T) gravity are possible. However, Harko et al. gave 
three classes of these models 

R+2f(T), 

f(R,T)={ /i(i2) + / 2 (n 

f 1 (R) + f 2 (R)f 3 (T). 

In this paper we are focussed to the first class, i.e. f(R,T) = R + 2f(T). 
For this model the field equations become 



R,», - ^Rg^u = kT^ v + 2f{T)T tMV + 



f(T) + 2pf'(T) 
where prime represents derivative with respect to T. 



9^, (9) 
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3 Exact Solutions of Bianchi Type / Universe 



In this section, we shall find exact solutions of Bianchi I spacetime in f(R, T) 
gravity. For this purpose, we use natural system of units (G = c = 1) and 
f(T) = XT, where A is an arbitrary constant. For Bianchi type I spacetime, 
the line element is given by 



ds 2 = dt 2 - A 2 (t)dx 2 - B 2 {t)dy 2 - C 2 {t)dz 2 , 



(10) 



where A, B and C are defined as cosmic scale factors. The Bianchi / Ricci 
scalar turns out to be 



R 



ABC 
A + B + C 



AB BC CA 
AB + BC + CA 



where dot denotes derivative with respect to t. 

Using Eq.()9]), we get four independent field equations, 



(11) 



AB 


BC 


CA 




+ 3X)p - Xp, 


AB + 


BC + 


CA 


= (8tt 


B 


C 


BC 


= Xp- 


- (8tt + 3A)p, 


B 


+ C + 


BC 


C 


A 


CA 


= Xp- 


- (8tt + 3A)p, 


c 


+ A + 


AC 


A 


B 


AB 




- (8tt + 3A)p. 


A 


+ B + 


AB 


= Xp 



(12) 
(13) 
(14) 
(15) 



These are four non-linear differential equations with five unknowns namely 
A, B, C, p and p. Subtracting Eq. (flij) from Eq.flTJ]), Eq.(Tl5D from Eq.dE 
and Eq.( fl5i) from Eq.( Tl2|) . we get respectively 



A B C(A B 
A~ B + C\A~ B 
B C A(B C 
B~C + A\B~C 
A C B f A C 
A~C + B \A ~ C 



0. 

o, 

: 0. 



(16) 
(17) 
(18) 
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These equations imply that 



B . 

— = d\ exp 
A 

— = d 2 exp 

— = d 3 exp 



dt 



dt 
c 2 I -s 



di 
cs / -j 



(19) 
(20) 
(21) 



where ci, c 2 , C3 and di, d 2 , cfe are integration constants which satisfy the 
following relation 

ci + c 2 + c 3 = 0, d 1 d 2 dz = l. (22) 

Using Eqs. fll9p - (l2Tl) . we can write the unknown metric functions in an explicit 

way 



dt 

9i / -5 



dt 
92 I 



where 
and 



A = api exp 
£? = aj? 2 exp 
C = aj?3 exp 

Pi = (dl~ 2 d3 -1 ) S , P2= (dl^fe -1 )*, P3 = (dlC?2 2 )5 



dt 
93 / 



ft 



2ci + c 2 ci — c 2 c x + 2c 2 

, 92 = — = — , Q3 = 



(23) 
(24) 
(25) 

(26) 
(27) 



3 ' " 3 ' ' 3 

It is mentioned here that p±, p 2 , £>3 and gi, g 2 , #3 also satisfy the relation 



P1P2P3 = 1, <?i + 92 + <?3 = 0. 



(28) 



3.1 Some Important Physical Parameters 

Now we present some important definitions of physical parameters. The 
average scale factor a and volume scale factor V are defined as 



a = VABC, V = a 3 = ABC. 



(29) 
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The generalized mean Hubble parameter H is given by 



H = -(H 1 + H 2 + H 3 ), 



(30) 



where Hi 



Ho 



{ . +i 2 — §, H 3 — ^ are defined as the directional Hubble 
parameters in the directions of x, y and z axis respectively. The mean 
anisotropy parameter A is 



.4 



1 3 



Hj~ H 
H 



(31) 



The expansion scalar 6 and shear scalar a 2 are defined as follows 



ABC 
A + B + C 1 



— <J uu (J 



where 



cr. 



.4 
A 

1 

: 2' 



AB BC CA 



AB BC CA 



-eh 



(32) 
,(33) 



(34) 



hfxu — 9/j.u ~ u^u v is the projection tensor. 



The deceleration parameter q is the measure of the cosmic accelerated 
expansion of the universe. It is defined as 



aa 

q = ~. 

cr 



(35) 



The behavior of the universe models is determined by the sign of q. The 
positive value of deceleration parameter suggests a decelerating model while 
the negative value indicates inflation. Since there are four equations and five 
unknowns, so we use a well-known relation [33] between the average scale 
factor a and average Hubble parameter H to solve the equations, 



H = la 



(36) 



where I and n are positive constants. 
Using Eqs.( l30|) and (|36|) . we get 



la 



l-n 



(37) 
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and the deceleration parameter becomes 

q = n-l. (38) 
Integrating Eq. (!37I) . it follows that 

a = {nlt + k 1 )^, n^O (39) 

and 

a = k 2 exp(lt), n = 0, (40) 

where k\ and k 2 are constants of integration. Thus we get two different 
models of the universe corresponding to these values of the average scale 
factor. 



3.2 Singular Model of the Universe 

Here we investigate the model of universe when n ^ 0, i.e., a = (nit + ki)n . 
In this case, the metric coefficients A, B and C takes the form 

qi(nlt + k\)~ 



A = p 1 (nlt + ki)™ exp 
B = p2{nlt + ki) » exp 
C = Ps(nlt + ki) » exp 



Z(n-3) 
l(n - 3) 



Z(n-3) 

The directional Hubble parameters Hi (i = 1,2,3) turn out to be 

Hi = 1 + Qi 

nlt + h (nlt + h)^' 

The mean generalized Hubble parameter and volume scale factor are 

/ 



(41) 
(42) 
(43) 

(44) 



H 



nit + k\ 

The mean anisotropy parameter becomes 



V = (nit + ki) 



A 



qi 2 + g2 2 + <?3 2 
3l 2 (nlt + k i y & - 2n y n ' 



(45) 



(46) 



The expansion scalar and shear scalar for this model are given by 

3/ 
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a 



2 <?i 2 + Q2 2 + qs 2 



nlt + ki " 2{nlt + A4) 6 /™ ' 
Using Eqs. f|T2|) - f|T5|) . the energy density of the universe is 



(47) 



12(A + 27r)(A + 47r) 



4(A + 3tt) 



3Z 2 



+ 



?l52 + <?2<?3 + <?3<7l 



3Z 2 (l-n) | qi z + q2 z + q ^ 



(nit + ki 



(nit + hi)' 



(4£ 



while the pressure of the universe becomes 



-1 



P 



12(A + 2vr)(A + 47r) 



47T 



3/ 2 



+ 



QlQ2 + <?2<?3 + <?3<?1 



(nlt + h) 2 (nlt + h) 1 . 



,, x s ,j 3/ 2 (l-n) gi 2 + g 2 2 + g 3 2 
+ (3A + 8tt)<; ^ ; - ^ + £ 



(nZt + A4) 2 



(49) 



3.3 Non-singular Model of the Universe 



For this model, n = and the average scale factor a = k20xp(lt) turns the 
metric coefficients A, B and C into 



A = p±k 2 exp(lt) exp 
£? = p 2 ^2 exp(Zt) exp 
C = exp(lt) exp 



<7i exp(— 3Zi) 

3/^ 
g 2 exp(-3/t) 

3/^ 
g 3 exp(-3/t) 



3lk 2 3 



The directional Hubble parameters Hi become 



Hi = I H ^exp(— 3/t). 

for 



(50) 
(51) 
(52) 

(53) 



The mean generalized Hubble parameter and volume scale factor turns out 
to be 

H = l, V = k 2 3 exp(3lt). (54) 
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The mean anisotropy parameter, expansion scalar and shear scalar are 



.4 



9i 2 + 92 2 + 93 2 
3/ 2 A; 2 6 exp(6^)' 



= 3L a 2 



qi 2 + g.2 + 93 2 

2k 2 6 exp(6lt) 



(55) 



For this model, the energy density and pressure of the universe takes the 
form 



P = 



V 



1 



12(A + 27r)(A + 47r) 

2j „2 

A<>3Z 2 +" 



A(X + 3n)\3l^ qiq2 + q2q3 + q3qi 



k 2 b exp(6lt) 



,2 , 91+92' 



Q3 



k 2 b exp(QU) 
-1 

12(A + 27r)(A + 4vr) 
+ (3A + 87r)|3/ 2 + 



4tt<^ 3l 2 + 



qiq2 + <?2<?3 + 9391 



fc 2 6 exp(6/t) 



2 , 9l92 + 9293 + 9391 



fc 2 b exp(6/t) 



(56) 



(57) 



4 Exact Bianchi Type V Solutions 

Here we shall explore Bianchi type V solutions in the context of f(R,T) 
gravity. The metric for the Bianchi type V spacetime is 



ds 2 



dt 2 - A 2 (t)dx 2 - e 2mx [B 2 (t)dy 2 + C 2 (t)dz 2 



(58) 



Here A, B and C are also cosmic scale factors and m is an any constant. 
The Ricci scalar for this spacetime is 



Using Eq.( 



R = -2 
D, we get 



A B 
A + B 



C 3m 2 AB BC CA 



+ 



+ 



+ 



C A 2 AB BC CA 



3m 2 



AB BC CA 
~AB + B~C + CA~^ 2 ~' 
B C BC _m 2 
B + C + ~BC~A 2 ~ Z 
C A CA _m 2 
C + A + AC~A 2 ~~~ 
A B AB m 2 
A + B + AB ~ I 2 " 



= (87T + 


3X)p - Xp, 


-Xp- 


(8vr + 3A)p, 


-Xp- 


(8tt + 3A)p, 


= Xp- 


(8tt + 3X)p 



(59) 

(60) 
(61) 
(62) 
(63) 
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and the 01-component turn out to be 



A B _C 



0. 



(64) 



We adopt the same procedure as for the Bianchi type I solutions to solve these 
equations. Here the equations Eqs.(|T6l)-(|T8l) are same as obtained previously 
but by making use of Eq. (1641 , we get the constraint equations as follows 



Pi = 1, p 2 =Ps = P, qi = 0, q 2 = -q 3 = Q- 
Thus, the metric coefficients become 



A = a, B = aP exp 



Ql d 4 



C = aP 1 exp 



-Ql d 4 



(65) 



(66) 



4.1 Singular Model of the Universe 

For the model of the universe when n^O, the metric functions A, B and C 
become 



A = (nit + 

B = P(nlt + kx)™ exp 



Q(nlt + ki) n 



C = p-\nlt + A;i)™exp 



Z(n-3) 
-Q(nlt + h)^ 



n^3 



l(n - 3) 

The directional Hubble parameters H\, H2 and H3 take the form 

/ 



#1 
H 2 



nit + ki 
I 



+ 



Q 



nlt + ki (nlt + k^n 
I Q 



nlt + kx (nlt + ki)n 
The mean anisotropy parameter becomes 

2Q 2 



(67) 
(68) 

(69) 



(70) 
(71) 

(72) 



.4 



3(nlt + A; 1 )( 6 - 2n )/ n ' 
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(73) 



The shear scalar for this model is given by 

Q 2 



a 2 = 



(nlt + h) 6 / n ' 



(74) 



It is to be noticed that the mean generalized Hubble parameter H, expansion 
scalar 9 and the volume scale factor V are same as in the case of Bianchi 
type / spacetime. The energy density and pressure of Bianchi V universe for 
this model turns out to be 



1 



- A 



12(A + 27r)(A + 47r) 
3l 2 (l-n) 



4(A + 3tt) 



2Q 2 



3/ 2 



Q 2 



(nlt + h) 2 (nlt + h)' 



(nlt + h) 2 



(nit + ki) n 



(75) 



-1 



P = 



12(A + 27r)(A + 47r) 
, 3l 2 (l -n 

+ (3A + 8tt) 



47T 



3l 2 



Q 2 



+ 



(nlt + h) 2 (nlt + h)^ 
2Q 2 



(nit + h) 2 (nit + k x ) 



(76) 



4.2 Non-singular Model of the Universe 

For the model when n = 0, the metric coefficients A, B and C turn out to 
be 



A = /c 2 exp(/t), 
B = Pk 2 exp(lt) exp 

C = P~ 1 k 2 exp(lt) exp 



Q exp(Slt) 

3lk? 
Q exp(Slt) 
3lh 3 



(77) 
(78) 

(79) 



The directional Hubble parameters Hi, H 2 and H 3 are 

Qexp(Slt) Qexp(Slt) 



Hi = l, H 2 = l + 



k 2 3 



H,= 



k 2 3 



(80) 



«2 K 2 

The mean anisotropy parameter and shear scalar for this model become 

2Q 2 2 Q 2 



A = 



3l 2 k 2 6 exp(6/t) 



k 2 6 exp(6lt) 



(81) 
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Here we also get the same volume scale factor V, expansion scalar 9 and 
mean generalized Hubble parameter H as shown in Eqs. fj53j) -f l54p . The energy 
density and pressure of universe here become 



P 



12(A 



27r)(A + 47r) 
A; 2 6 exp(6/t) JJ' 



4(A + 3ttW 3r - 



Q 2 



k 2 b exp(6lt) 



(82) 



V 



12(A + 27r)(A + 47r) 
+ (3A + 87r)i 3Z 2 + 



4ir< 31 



Q 2 



k 2 exp(6lt) 



2Q 2 



k 2 b exp(Qlt) 



(83) 



5 Concluding Remarks 



This paper is devoted to discuss the current phenomenon of accelerated ex- 
pansion of universe in the framework of newly proposed f(R,T) theory of 
gravity. For this purpose, we take f(R, T) = R + 2 AT and explore the exact 
solutions of Bianchi types I and V cosmological models. We obtain two ex- 
act solutions for both spacetimes using the assumption of constant value of 
deceleration parameter and the law of variation of Hubble parameter. The 
obtained solutions correspond to two different models of universe. The first 
solution forms a singular model with power law expansion while the second 
solution gives a non-singular model with exponential expansion of universe. 
The physical parameters for both of these models are discussed below. 

The singular model of the universe corresponds to n / with average 
scale factor a = (nit + fci)™. This model possesses a point singularity when 
t = t s = — ^j. The volume scale factor vanishes and the metric coefficients 
A, B and C vanish at this singularity point. The cosmological parameters 
H\, H 2 , Hs, H, 9, and a 2 are all infinite at this point of singularity. The 
mean anisotropy parameter A also becomes infinite at this point for < n < 3 
and vanishes for n > 3. The energy density and pressure of universe are also 

2 

infinite at this epoch. Moreover, the isotropy condition, i.e., ^- — > as 
t — > oo, is verified for this model. All these conclusive observations suggest 
that the universe starts its expansion with zero volume, strong pressure and 
energy density from t = t s and it will continue to expand for < n < 3. 
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Now we discuss the non-singular model of the universe corresponds to 
n — 0. For this model the average scale factor is a = k 2 e^{lt). The non- 
singularity is due to the exponential behavior of the model. The expansion 
scalar 6 and mean generalized Hubble parameter H are constant in this 
case. For finite values of t, the physical parameters H-y, H 2 , H 3 , a 2 and 
A are all finite. The metric functions are defined for finite time and the 
isotropy condition is satisfied. The pressure and energy density of universe 
become infinite in the limiting case when t — > — oo. This shows that the 
universe evolved from an infinite past with a massive energy density and 
pressure. There is an exponential increase in the volume as the time grows. 
This shows that the universe started its expansion a long time ago with zero 
volume. 
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